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ABSTRACT 



We report the detection in emission of various highly excited rotational transitions of HNC (J = 6-5 through J =12-1 1) toward the 
carbon-star envelope IRC +10216 using the HIFI instrument on-board the Herschel Space Observatory. Observations of the J = 
1-0 and J = 3-2 lines of HNC with the IRAM 30-m telescope are also presented. The lines observed with HIFI have upper level 
energies corresponding to temperatures between 90 and 340 degrees Kelvin, and trace a warm and smaller circumstellar region than 
that seen in the interferometric maps of the 7=1-0 transition, whose emission extends up to a radius of 20". After a detailed chemical 
and radiative transfer modeling, we find that the presence of HNC in the circumstellar envelope of IRC +10 216 is consistent with 
formation from the precursor ion HCNH + , which in turn is produced through several proton transfer reactions which are triggered by 
the cosmic-ray ionization. We also find that the radiative pumping through A 2 1 yum photons to the first excited state of the bending 
mode v 2 plays a crucial role to populate the high-./ HNC levels involved in the transitions observed with HIFI. Emission in these 
high-./ rotational transitions of HNC is expected to be strong in regions which are warm and dense and/or have an intense infrared 
flux at wavelengths around 21 /im. 
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1. Introduction 

Hydrogen isocyanide (HNC), a metastable isomer of HCN ly- 
ing 0.6 eV higher in energy, is ubiquitous in the interstellar 
medium. First observed in several objects through a line at 90.7 
GHz (e.g. Zuckerman et al. 1972| l, the definitive confirmation 
of this line, the J = 1-0 rotational transition of HNC, came 
through spectroscopic laboratory experiments in the late 1970s 
(Blackman et al. 1976). HNC is widely observed in different 
types of astronomical regions, such as diffuse clouds fLiszt &| 
|Lucas 2001)1, cold dark clouds ( Hirota et al. 1998 1, star-forming 



regions (Schilke et al. 



evolved stars dBujarraba 



The HIFI instrument on-board the Herschel Space 
Observatory, operating in the 480-1910 GHz frequency 
range, allows to cover high-energy rotational transitions of 
HNC. In this study, we present the detection in the C-star 
envelope of IRC +10216 of the J = 6-5 through J = 12-11 
rotational transitions of HNC. Observations of the 7=1-0 and 
J = 3-2 transitions obtained with the IRAM 30-m telescope 
are also presented. The large set of observed lines allows us to 
investigate in detail the abundance distribution and excitation of 
HNC in IRC +10216's circumstellar envelope. 



1992), circumstellar envelopes around 



et al.| |l994f, exter nal galaxies (|A~alto et| 2. Observations and results 



[aT 2002), and even in high redshift objects (Guelin et al. 2007 ). 
The main formation mechanism of HNC is thought to be the 
dissociative recombination of the ion HCNH + ([Semaniak et al. 
|2001| l. HNC has mostly been observed only in the J = 1-0 tran- 
sition at ~3.3 mm, preventing to carry out a detailed excitation 
analysis (Sarras in et al. 2010| >. 



* Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 

** Based on observations carried out with the IRAM 30-meter tele- 
scope. IRAM is supported by INSU/CNRS (France), MPG (Germany) 
and IGN (Spain). 



The HIFI observations were obtained in the context of a high res- 
olution line survey of IRC +10216 ( Cernicharo et al. [2010). The 
data were taken in double beam-switching mode with a spectral 
resolution of 1.1 MHz and channel spacings of 0.5 MHz, and 
processed using the standard Herschel pipeline up to Level 2, 
providing fully calibrated spectra. The HNC lines where all ob- 
served around the maximum luminosity, i.e. phase = 0. For 



details regarding the reduction procedure we refer to Cernicharo 
|et al.| ( |2010| l. The total integration time per frequency setting 
ranged from 2 to 14 minutes, resulting in antenna temperature 
rms noise levels of 8-15 mK per 1.1 MHz channel for the HNC 
spectra in bands la, 2a, 2b, and 3a, and of 50 mK for the spectra 
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Table 1. HNC observed line parameters 
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Fig. 1. HNC rotational lines observed toward IRC +10216 with 
the IRAM 30-m telescope (7 = 1-0 and 7 = 3-2) and with HIFI 
(7 — 6 - 5 to 7 — 12-11). The red and blue curves correspond 
to the models described in Sect.[3]and[4] 



above 900 GHz. At the frequencies of the HNC lines observed 
by HIFI (7 = 6-5 through 7 = 12-11), the beam sizes range from 
39" down to 20", and is therefore comparable to the size of the 
emitting region of HNC in IRC +10216. The 7 = 1-0 transi- 
tion, as mapped with the IRAM Plateau de Bure Interferometer 
(PdBI), shows a ring-like distribution extending up to a radius 
of nearly 20"(see Guel in et al.||1997| >, and the higher-/ transi- 
tions observed with HIFI are expected to have a more compact 
distribution. 

The 7=1-0 transition of HNC was observed before 2009 
with the IRAM 30-m telescope in the course of a A 3 mm spec- 
tral line survey of IRC +10216 (Cernicharo et al., in prep.). The 
achieved T* A rms noise level after averaging different spectra at 
the frequency of the 7=1-0 line is 3 mK per 1 MHz chan- 
nel. The 7 = 3-2 line was observed in January 2010 during a 
A 0.9 mm spectral line survey of IRC +10216 (Kahane et al., 
in prep.), using the new EMIR receivers. The line was observed 
with a spectral resolution of 2 MHz during about 3 h of good at- 
mospheric conditions (7 sys = 800 K) resulting in a FJ rms noise 
level of 6 mK per channel. Additionnaly, we make use of the 
PdBI observations of the 7=1-0 line. These observations are 
described in |Guelin et al.| l fl997) and references therein. 
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The observed lines of HNC are shown in Fig. [T] and the line 
and telescope parameters are given in Table [T] All the lines ap- 
pear free of contamination from other lines. The 7=1-0 and 
7 = 3-2 lines observed with the IRAM 30-m telescope are the 
most intense ones and show a U-shaped profile, typical of opti- 
cally thin emission resolved by the telescope. On the other hand, 
the lines observed with HIFI show a flat-topped profile, which is 
typical of optically thin lines not resolved by the telescope, i.e. 
geometrically diluted in the main beam of HIFI. 

3. Modeling 

To learn about the excitation conditions and the formation mech- 
anism of HNC in the envelope of IRC +10216 we have carried 
out modeling tasks. The model, basically taken from Fonfria et 
|aL| ( |2008) and Agundez ( 2009 [ l, considers the presence of a cen- 
tral star with an effective temperature T* = 2330 K and a radius 
R t = 4 x 10 13 cm, surrounded by a spherical envelope of gas 
and dust. The adopted mass loss rate is 2 x 10~ 5 Mo yr _1 , with 
a distance of 130 p c, and a terminal expansio n velocity of 14.5 
km s . Following Cordiner & Millar (20091, we consider the 
presence of density-enhanced shells. 

The gas temperature profile was derived using 12 CO and 
13 CO HIFI and IRAM 30-m observations. The adopted 12 CO 
abundance relative to H2 was 6 x 10~ 4 and the 12 CO/ 13 CO abun- 
dance ratio 45 (Kah ane et al.||2000l ). The best fit to the whole set 
of lines is given by a 3-fold power law, i.e. 7* oc (r/R m )~" with 
n = 0.55 for r <75R t ,n = 0.85 for r < 200 fl, and n = 1.4 for 
greater radii, with the scale fixed by setting 7* = T* at r - R„. 
The estimated uncertainty in 7* is of 50 %. The physical param- 
eters that describe the circumstellar envelope, i.e. the dust and 
gas temperature as well as the H2 volume density are presented 
in Fig. 

As will be shown below, infrared (IR) pumping to excited vi- 
brational states of HNC plays an important role in the excitation 
of the observed lines. Hence, a good description of the IR flux is 
needed. For that purpose the spectral energy distribution (SED) 
observed by ISO at phase <p = 0.24 (Cernicharo et al. 1999[ ) has 
been modeled to derive dust opacities suitable for IRC +10216, 
which are given by a mixture of silicon carbide, amorphous car- 
bon and MgS dust grains (see Fig. |3j. The HNC lines were ob- 
served at phase <p = 0.95-0.05, around the maximum luminosity. 
Therefore, the dust temperature and SED were recomputed mod- 
ifying the star radius accordingly and using photometric points 
observed at maximum ( |Le Bertre) 1 1 992| l . More details concern- 
ing the modeling of IRC +10216's SED can be found in |De 
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Fig. 2. Physical parameters (i.e. dust and gas temperatures, H2 
volume density) used in the present work to describe the cir- 
cumstellar envelope, density-enhanced shells are introduced in 
the current modeling as can be seen in the H2 volume density 
curve. The first shell is located at 15" and we assume an inter- 
shell distance of 12", the shells being 2" wide. 
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Fig. 3. Comparison between the SED of IRC +10216 as mod- 
eled and as observed by ISO at <p — 0.24. Wavelengths of some 
vibrational transitions of HNC are indicated by vertical lines. 
The vibrational transition indicated at ~ 11 /urn correspond to 
the 000 — > 020 vibrational band and the transition at ~ 2 1 /mi to 
the 000 -> 010 and 010 -> 020 bands. The dotted vertical line 
around 45 yum indicates the ISO LWS / SWS division. The inset 
panel shows the computed SED at different phases. 



Beck et aT7| ( 2012| l. The parameters used in the present study to 
describe the dust (i.e. the dust opacities and dust-temperature 
profile) are the same as those described in the latter article. The 
differences in the computed SEDs mainly come from the inclu- 
sion of the density-enhanced shells (see Sect. |4j, of which the 
main effect is that the flux at wavelengths longer than 25 um is 
slightly increased. 

3.1. Chemical model 

The chemical composition of the gas has been computed as it ex- 
pands from a radius of 2 x 10 14 cm up to 10 18 cm. To deal with 
the density-enhanced shells we have proceeded as Cordi ner"&| 
Millar ( 2009), running a model for the smooth component of the 
envelope and another for a density-enhanced shell. The cosmic 
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Fig. 4. Abundance distribution of HNC obtained with the chemi- 
cal model (red curve). The modification introduced to reproduce 
the observations (blue lines in FigJT] ) is shown by to the blue 
curve. 
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Fig. 5. Abundances of some of the species (i.e. HCNH + , HCN, 
NH3 and e-) related to the HNC formation as obtained from the 
standard chemistry model. 



|Cernicharo [20 06) and the ultraviolet interstellar radiation field 
of Draine| ( 1978| l is adopted. The chemical network comprises 
470 gas phase species (composed of H, He, C, N, O, Si, P, and 
S) linked by 7400 reactions, whose rate constants are taken from 
databases such as udfaO^t Woodall et al. 2007 1 and from the 
recent literature on chemical kinetics. The abundances at the ini- 
tial radius are taken from Agundez (2009]!. 

The resulting abundance distribution of HNC in the enve- 
lope of IRC +10216 is plotted in Fig. [4] and the abundance of 
the species chemically related to HNC in Fig. [5] Apart from the 
abrupt changes at 15", 29", and 43", which correspond to the 
presence of density-enhanced shells, we see that the fractional 
abundance of HNC starts to increase at small radii, well before 
the molecular shell located at 15". This occurs because the for- 
mation of HNC is, to a large extent, driven by cosmic rays that 
can penetrate deeper in the envelope than interstellar ultravio- 
let photons. In the regions inward of 10 16 cm, HNC is mostly 
formed through the reaction: 



HCNH + + NH 3 



hnc + nh; 



(1) 



Ammonia is relatively abundant in the inner layers o f 
IRC +10216 (2 x 10~ 6 relative to H 2 ; |Hasegawa et al.||2006] >, 



ray ionization rate of H2 is taken as 1.2 x 10 s (Agundez & 



See http://www.udfa.net 
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while the precursor ion HCNH + is formed from several reac- 
tions of proton transfer to HCN (a species which is also fairly 
abundant in the inner envelope of IRC +10216, 2 x 10" 5 relative 
to H2 according to |Fonfrfa et al.|20 08 ): 

HCO + + HCN -> HCNH + + CO (2) 



where HNC has a significant abundance, the collisional excita- 
tion of the vibrational levels of HNC is expected to be negligible 
(c.f. Sect. [4}. Therefore, for ro-vibrational transitions, the same 
set of rate coefficients was used but they were scaled down by an 
arbitrary factor of 100. 



C 2 H+ + HCN -> HCNH + + C 2 H 2 



( 3 ) 4. Results 



Hj + HCN 



HCNH + + H 2 



(4) 



These processes are driven by the cosmic-ray ionization of H 2 
and the subsequent fast formation of H^, which yields HCO + 
and C 2 H^ through proton transfer to CO and C 2 H 2 , respectively. 
Reactions ([TJ-^jl have been experimentally found to be rapid (rat e 



constants are in excess of 10 9 cm J s 1 at 300 K; Anicich|2003 1 



For reaction ([TJ, however, it is not clear what are the branching 
ratios yielding HCN and HNC, so that we consider that both 
channels occur at the same rate, as is done by e.g. the databases 
udfa06' and KIDApl This assumption will be discussed in Sect. 

The peak abundance of HNC, 7 x 10 8 relative to H 2 , is 
reached beyond 10 16 cm, concretely at about 4 x 10 16 cm. In 
this region the formation of HNC occurs mostly through the re- 
actions: 



HCNH + + e" -> HNC + H 



C + NH 2 -> HNC + H 



(5) 
(6) 



where the formation of HCNH + is again driven by cosmic rays 
while reaction |6]i occurs as a consequence of the photochem- 
istry, with atomic carbon and NH 2 coming from the photodisso- 
ciation of acetylene and ammonia, respectively. Reac tion (|5 ]l has 
been well studied in the laboratory (Seman iak et al.||2001) , but 
there are still some uncertainties regarding the branching ratios 
yielding HCN and HNC. On the other hand, the kinetics of re- 
action |6]) is poorly known and the rate constant has been taken 
from an estimate of Smith et al. (2004). The destruction of HNC 
in IRC +10216's envelope is dominated by photodissociation, 
for which the rate is taken from udfa06'. Despite the lack of 
accurate branching ratios and rate constants for some of the re- 
actions involved in HNC formation, its predicted abundance has 
the right order of magnitude, although these inaccuracies may 
introduce some errors in the abundance radial distribution calcu- 
lated, as will be shown in Sect. [4] 

3.2. Radiative transfer model 

The radiative transfer calculations were performed with the code 
described in |Daniel & Cernicharo| ( |2008| ). The first 20 rotational 
levels of HNC in the ground and in the excited vibrational states 
v 2 = 1,2 (bending mode at 21 fim with {-type doubling), vj = 1 
(NH stretching mode at 2.7 fim), and V3 = 1 (NC stretching 
mode at 4.9 yum) were included with the spectroscopic parame- 
ters given by Maki & Mellau (200T|. The experimental electric 
dipole moment of HNC, 3.05 Debye ( |Blackman et al.[ |1976) , 
was used for the rotational transitions in each vibrational state. 
For ro-vibrational transitions the dipole moments were taken 
from |Harris et al. ( 2002|l. The collisional rate coefficients com- 
puted by Dumouch el et al.| ( [2010) > for the first 26 rotational levels 
and for temperatures between 5 and 500 K were used for pure 
rotational transitions. For the densities prevailing in the region 



See http : //kida . obs . u-bordeauxl . f r 



The HNC line profiles that result from the abundance distribu- 
tion given by the chemical model are plotted in Fig. [T[ as red 
lines. It is seen that the model overestimates the intensity of 
some of the HIFI lines, in particular the J = 6 - 5. A better 
overall agreement is obtained if the abundance of HNC is low- 
ered to 7 x 10~ 9 , relative to H 2 , for r < 15" (see blue lines in 
Fig. [4] and Fig. [TJ. It, thus, seems that the chemical model over- 
estimates the abundance of HNC in the regions inward of 15" by 
a factor of ~ 5. 

4.1. HNC excitation 

The influence of the different vibrational excited states on the 
excitation of the observed HNC lines is illustrated in Fig. [6] It is 
seen that the line intensities are largely increased when the v 2 =l 
state is included, and to a lesser extent when the the stretching 
modes are considered. From the results shown in this figure, it is 
concluded that the excitation of HNC rotational levels is mainly 
dominated by infrared pumping to the excited vibrational states, 
and to a lesser extent by inelastic collisions. Therefore, results 
should not be substantially altered by the uncertainty in the 7* 
profile, something that has been checked adopting alternative 7* 
profiles. We note that this also validates the use of an arbitrary 
set of ro-vibrational collisional rate coefficients. This point has 
been checked by considering scaling factors that differ by a fac- 
tor 2 with respect to the value of 100 used to scale down the 
ro-vibrational rate coefficients. No differences were found be- 
tween the various results. Owing to the efficiency of the radiative 
pumping to the vibrational levels, the accuracy of the radiative- 
transfer model largely relies upon the correctness of the SED. In 
order to test the error that would result from inaccuracies in the 
SED modeling, we adopted alternate models for the dust com- 
position, that result in qualitatively good representations of the 
overall SED, but with variations in its shape (i.e. with maximum 
variations of the order of 10-20% at a given frequency, in com- 
parison to the current model). From these tests, it is concluded 
that the line fluxes will show variations of the same order as the 
variations observed for the continuum at 11 and 21 /im. Thus, 
even if the quality of the SED can play a role in explaining the 
differences observed between the observations and the model, 
the effect should be modest and may not explain the differences 
observed (a factor ~ 2 for the 7 = 6-5 line). 

4.2. Density-enhanced shells 

Another aspect worth considering is the inclusion of density en- 
hanced shells. Such shells have been characterized through opti- 
cal imaging by Mauron & Huggins ( 2000 ) and would originate 
from an episodic increase of the mass loss rate of the star. From 
this study, it appears that their exact characterization is complex: 
rather than continuous spherical shells, various disrupted arcs of 
matter are detected. These shells have more recently been char- 
acterized with the PACS instrument on-board Herschel ( |Decin| 
|et al. 201 1) . This study shows that the shell structure extend 
at least out to ~320" and that the shells have a density which 
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Fig. 7. Images of the HNC(7 =1-0) emission integrated over the 10 km/s central velocity around the Vlsr of the source. The 
fluxes are given in Jy/beam, with a beam size which correpond to the Plateau de Bure interferometer (~4"). The various panels 
correspond to (a) IRAM PdBI observations reported in Guelin et al. ( 1997] ) (b) model with density-enhanced shells and with a 
central decrease of the abundance (blue curve in Fig. |4ji (c) model with density-enhanced shells and with the abundance given by 
the standard chemical model (red curve in Fig. |4j (d) model without density-enhanced shells and with the abundance given by the 
standard chemical model (blue curve in Pig . [8j> . In panels (b), (c) and (d), the isocontours correspond to the IRAM PdBI observations 
displayed in panel (a). 



is typically enhanced by a factor 4 with respect to the inter- 
shell medium. An attempt in introducing such shells in order to 
deal with the analysis of molecular emission has been proposed 
by Cordiner & Millar (2009). These authors introduce 2" wide 
shells with an intershell distance of 12". This model is successful 
in reproducing the spatial distribution of C,,H molecular species, 



and we therefore adopted the same structure for the shells in the 
current work. 

The main effect introduced by these shells is on the ra- 
dial profile of the molecular abundance predicted by the chemi- 
cal modeling. Indeed, without density-enhanced shells, the pho- 
todissociation of the molecules is efficient down to smaller radii, 
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Fig. 6. Effect induced by the inclusion of the vibrational modes 
of HNC on the lines observed by the IRAM and Herschel tele- 
scopes, for the model with unaltered X(HNC) profile. 



as compared to a model including such shells, and therefore the 
peak abundance of the molecules formed by photochemistry is 
shifted inward. This result is similar to the one described by [Del 
Beck et al. ( 20 1 2 1 while dealing with the interpretation of the 
C2H HIFI observations. This is illustrated in Fig. [8]for the case 
of HNC where the abundance of the molecule is presented for 
the two models. The line profiles obtained for the two HNC ra- 
dial distributions are shown in Fig. [9] Except for the 7=1-0 
line, omitting the density-enhanced shells would result in a bet- 
ter agreement between the model and observations. Still, in the 
model without density-enhanced shells, the 7 = 6-5 line is 
overestimated, the intensity being ~ 35% higher in the model 
in comparison to the observation. Fig. [7] shows the flux of the 
7=1-0 line, integrated over a 10 km/s wide interval centered 
on the Vlsr of the source. Panel (a) corresponds to the IRAM 
Plateau de Bure observations described in Guel in et al.| ( |1997] l. 
In this figure, the extent of the emission predicted by the models 
is shown for various cases: (b) with density-enhanced shells and 
including a central drop in the HNC abundance, as derived from 
the HIFI observations (c) with density-enhanced shells, with the 
abundance profile predicted by the chemical modeling (d) with- 
out density-enhanced shells and with the abundance profile pre- 
dicted by the chemical modeling. First, a comparison of pan- 
els (c) and (d) shows that the inclusion of the density-enhanced 
shells has a profound influence on the extent of the emission. 
As commented previously, the density-enhanced shells attenu- 
ate the UV field, reducing the amount of photodissociation. This 
entails that the inclusion of the density-enhanced shells increases 
the extent of the emission of the 7=1-0 line to larger radii. 



NO enhanced shells 
enhanced shells 




20 40 
radius (arcsec) 



Fig. 8. HNC abundance as a function of radius, without (blue) 
and with (red) density enhanced shells. 



Comparing the prediction of the two models with the observa- 
tions, it seems that a good fit to the observations would corre- 
spond to a structure intermediate, in the sense of the attenuation 
of the external UV field, to the models with and without density- 
enhanced shells considered here. Moreover, none of these mod- 
els is able to reproduce the extent of the ring-shaped structure 
which is observed. On the other hand, and considering panel (b), 
we see that the central drop introduced to fit the HNC lines ob- 
served with HIFI permits to reproduce this ring-shaped struc- 
ture. Moreover, we note that such a decrease in the HNC abun- 
dance would have to be introduced in any model intermediate 
to the two models presented in Fig. [8] Indeed, the model with- 
out density-enhanced shells corresponds to the lowest abundance 
possible for HNC, as predicted by the chemical modeling, for 
radii r < 12" . Whatever could be the structure considered for the 
shells, the abundance of HNC would have to be reduced in this 
region in order to reproduce the observed 7 = 6-5 line and to fit 
the ring-shaped structure of the 7=1-0 emission. Finally, by 
considering various studies that deal with the geometrical struc- 
ture of IRC +10 216, it is evident that density-enhanced shells 
are present in the source and have to be in cluded (e.g.|Mauron &| 
|Huggins||2000l|Cordiner & Millar| |2009| |De Beck et al.||2012) . 
However, the model adopted in the present study is a simplifi- 
cation of the exact structure of the shell morphology which is 



observed through optical imagin g (|Mauron & Huggins] [2000 ) 
or through CO lines fFong et ah] |2003||Cernicharo et al.j |2012 1 



and this introduces a source of uncertainty in the current mod- 
eling. But, as discussed earlier when considering the "extreme 
case" of an envelope without density-enhanced shells, the fact 
that the HNC abundance has to be decreased with respect to the 
prediction of the chemical modeling should not be altered by the 
description adopted to model the shells. 

4.3. HNC chemistry 

As discussed in Sect. |3.1| the formation of HNC in the region 
r < 30" is mainly driven by three reactions, of which two have 
poorly constrained branching ratios or rate constants. In what 
follows, we examine how the uncertainties on these reactions, 
i.e. reactions (1) and (6) of Sect, 
chemical modeling overestimates t 



3.1 



could explain why the 
leTHNC abundance. Fig. [10] 
shows the HNC abundance profile obtained by modifying the 
rate constant of reactions (1) and (6). For the case of reaction 
(1), the alteration consists in modifying the branching ratio of 
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Fig. 9. HNC line profile for the model with density-enhanced 
shells (red) and for the model without shells (blue). 



the reaction and for the curve shown in Fig. 10 the abundance is 
shown for the case where 75% of the collisions lead to HCN and 
25% to HNC. For the case of reaction (6), the abundance shown 
in Fig. [T0| corresponds to the case where this reaction is omit- 
ted from the chemical network. For these two cases, the corre- 



sponding line profiles are presented in Fig. 1 1 First, we see from 
Fig. [lOjthat altering reaction (1) can introduce a reduction of the 
HNC abundance for radii below 15". This is the region where 
the HNC abundance has been found to be overestimated when 
modeling the observations. However, considering the line pro- 
files shown in Fig. 1 1 solely changing the branching ratio of this 
reaction would not enable to reproduce the observations. Indeed, 
fitting the 7 = 6-5 line would require to lower the amount of 
HNC produced via this reaction below 25% and in such a case, 
the higher excited lines would have their intensities largely un- 
derestimated. The removal of reaction (6) from the chemical net- 
work modifies the HNC abundance for radii 15" < r < 40". The 
effect on the line profile, in comparison to the standard network 
is rather low, and only lead to a slight lowering of the intensities 
of the 7=1-0 and 7 = 3-2 lines. So within to the current 
chemical network the most probable source of uncertainty lies in 
the chemical kinetics of reaction (1) and its branching ratio. We 
stress that the amount of HNC produced via this reaction also 
depends on the NH3 abundance. In the case of IRC +10216, the 
NH3 abundance in the innermost part of the envelop e has been 
derived to be in the range 10~ 7 <X(NH 3 ) < 2 x 10" 6 ( Hasegawa 
et al. 2006; Keady & Ridgway 1993] >. In the current model, we 
assume an abundance that corresponds to the upper limit (see 
Fig. [2]). We tested the lower limit for the NH3 central abundance 
and found that the decrease in the HNC abundance was lower 
than in the model where the branching ratio of reaction (1) is 
modified. 



4.4. Conclusion 

Several conclusions may be obtained from the radiative trans- 
fer modeling. Firstly, HNC transitions with higher energies have 
smaller emitting regions. For instance, the emission region of 
the high-7 transitions observed with HIFI is noticeably smaller 
than that of the 7=1-0 transition, which emits up to a radius 
of 20" (see Guelin et al. 1997). The model clearly illustrates 
this on Fig. 12 where the velocity-integrated intensity of vari- 
ous lines is plotted as a function of radius. A second conclusion 
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Fig. 10. HNC abundance profile derived for various cases. The 
red curve correspond to the abundance obtained with the stan- 
dard chemical network. The blue curve is obtained by changing 
the branching ratio of reaction (1) of Sect. 3.1 and the dashed 



blue curve removing reaction (6) from the chemical network (see 
text). 
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Fig. 11. HNC line predicted for the abundance profiles shown 
in Fig. 10 The red and blue lines correspond respectively to an 



alteration of reaction (1) or (6). 

is that infrared pumping to the vi=l state is the main mechanism 
to populate the HNC rotational levels of the ground vibrational 
state in IRC +10216's envelope. Thus, the emission of high-7 
lines of HNC is favored in regions which are warm and/or have 
a large radiation flux at A 21 /zm. It is worth here to comment 
on the detection of HNC V2=l in the carbon-rich pre-planetary 
nebula CRL 618 ( |Schilke et aL]|2003] ). These authors found that 
vibrationally excited HNC is confined to the hot circumstellar 
layers, although it is likely pumped by A 21 //m photons be- 
cause of the high radiation flux at this wavelength in CRL 618. 
In IRC +10216, the radiative transfer model predicts that ro- 
tational lines within the V2=l state would have intensities of a 
few mK, and thus may have escaped detection due to the lack of 
sensitivity. 



5. Summary 

The detection of highly excited rotational lines of HNC toward 
IRC +10216 evidences the existence of a warm HNC com- 
ponent whose emission arises from regions smaller (extending 
up to a few arcseconds) than those traced by the fundamen- 
tal rotational transition, previously mapped at A 3 mm with the 
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Fig. 12. Velocity integrated intensity as a function of impact pa- 
rameter for various HNC lines. 



PdBI, and whose emission extends up to 20" (see |Guelin et al. 



1997 1. These results are still compatible with HNC being formed 
not close to the star, but in the intermediate and outer layers 
where cosmic-ray ionization drives the formation of the precur- 
sor ion HCNH + . This study has shown that the infrared pumping 
through the first excited state of the bending mode V2 dominates 
the excitation of the high-/ rotational levels involved in the tran- 
sitions observed with HIFI. Emission in these high-/ rotational 
transitions of HNC is expected to be strong in regions where 
HNC is either present in a warm and dense region and/or sur- 
rounded by intense infrared radiation at wavelengths around 21 
fim. 

The present study shows that the HNC abundance can be 
qualitatively reproduced by the chemical modeling and using 
the currently available reaction network. However, a quantitative 
estimate of the HNC abundance points toward multiple origins 
concerning the discrepancies between observations and the mod- 
els. The discrepancies may arise from a too simplistic descrip- 
tion of the physical structure of the envelope and particularly 
of the density-enhanced shells. Another origin could be in the 
description of the SED of the source or, finally, in the uncertain- 
ties in the reaction rates, branching ratio and missing reactions 
in the chemical network. We stress that a 2D description of the 
chemistry and radiative transfer modeling could enable to disen- 
tangle between these various possibilities but such a work would 
require additional observations in order to constrain the geomet- 
rical structure of the source. In the next years, the observing ca- 
pabilities that will be available through the ALMA interferom- 
eter will shed light on this, since the J=l-0, 3-2, 4-3, 7-6 and 
9-8 HNC lines could then be observed in this object with high 
angular resolution. 
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